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FIELD OF THE INVENTION 
[0 002] The present invention generally relates to a method of estimating electncal 
parameters of an earth formation with a simplified measurement device model. 
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BACKGROUND OF THE INVENTION 
[0003] Typical petroleum drilling operations employ a number of techniques to gather 
information about earth formations during and in conjunction with drilling operations such as 
Wireline Logging, Measurement- While-Drilling (MWD) and Logging- While-Drilling (LWD) 
operations. Physical values such as the electrical conductivity and the dielectric constant of 
an earth formation can indicate either the presence or absence of oil-bearing structures near a 
drill hole, or "borehole." A wealth of other information that is useful for oil well drilling and 
production is frequently derived from such measurements. Originally, a drill pipe and a drill 
bit were pulled from the borehole and then instruments were inserted into the hole in order to 
collect information about down hole conditions. This technique, or "wireline logging," can 
be expensive in terms of both money and time. In addition, wireline data may be of poor 
quality and difficult to interpret due to deterioration of the region near the borehole after 
drilling. These factors lead to the development of Logging- While-Drilling (LWD). LWD 
operations involve collecting the same type of information as wireline logging without the 
need to pull the drilling apparatus from the borehole. Since the data are taken while drilling, 
the measurements are often more representative of virgin formation conditions because the 
near-borehole region often deteriorates over time after the well is drilled. For example, the 
drilling fluid often penetrates or invades the rock over time, making it more difficult to 
determine whether the fluids observed within the rock are naturally occurring or drilling 
induced. Data acquired while drilling are often used to aid the drilling process. For example, 
MWD/LWD data can help a driller navigate the well so that the borehole is ideally positioned 
within an oil bearing structure. The distinction between LWD and MWD is not always 
obvious, but MWD usually refers to measurements taken for the purpose of drilling the well 
(such as navigation) whereas LWD is principally for the purpose of estimating the fluid 
production from the earth formation. These terms will hereafter be used synonymously and 
referred to collectively as "MWD/LWD." 

[0004] In wireline logging, wireline induction measurements are commonly used to gather 
information used to calculate the electrical conductivity, or its inverse resistivity. See for 
example U.S. Patent 5,157,605. A dielectric wireline tool is used to determine the dielectric 
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constant and/or resistivity of an earth formation. This is typically done using measurements 
which are sensitive to the volume near the borehole wall. See for example U.S. Patent 
3,944,910. In MWD/LWD, a MWD/LWD resistivity tool is typically employed. Such 
devices are often called "propagation resistivity or "wave resistivity" tools, and they operate 
at frequencies high enough that the measurement is sensitive to the dielectric constant under 
conditions of either high resistivity or a large dielectric constant. See for example U.S. 
Patents 4,899,1 12 and 4,968,940. In MWD applications, resistivity measurements may be 
used for the purpose of evaluating the position of the borehole with respect to boundaries of 
the reservoir such as with respect to a nearby shale bed. The same resistivity tools used for 
LWD may also used for MWD; but, in LWD, other formation evaluation measurements 
including density and porosity are typically employed. 

[0005] For purposes of this disclosure, the terms "resistivity" and "conductivity" will be 
used interchangeably with the understanding that they are inverses of each other and the 
measurement of either can be converted into the other by means of simple mathematical 
calculations. The terms "depth," "point(s) along the borehole," and "distance along the 
borehole axis" will also be used interchangeably. Since the borehole axis may be tilted with 
respect to the vertical, it is sometimes necessary to distinguish between the vertical depth and 
distance along the borehole axis. Should the vertical depth be referred to, it will be explicitly 
referred to as the "vertical depth." 

[0006] Typically, the electrical conductivity of an earth formation is not measured directly. 
It is instead inferred from other measurements either taken during (MWD/LWD) or after 
(Wireline Logging) the drilling operation. In typical embodiments of MWD/LWD resistivity 
devices, the direct measurements are the magnitude and the phase shift of a transmitted 
electrical signal traveling past a receiver array. See for example U.S. Patents 4,899,112, 
4,968,940, or 5,811,973. In commonly practiced embodiments, the transmitter emits 
electrical signals of frequencies typically between four hundred thousand and two million 
cycles per second (0.4 - 2.0 MHz). Two induction coils spaced along the axis of the drill 
collar having magnetic moments substantially parallel to the axis of the drill collar typically 
comprise the receiver array. The transmitter is typically an induction coil spaced along the 
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axis of a drill collar from the receiver with its magnetic moment substantially parallel to the 
axis of the drill collar. A frequently used mode of operation is to energize the transmitter for 
a long enough time to result in the signal being essentially a continuous wave (only a fraction 
of a second is needed at typical frequencies of operation). The magnitude and phase of the 
signal at one receiving coil is recorded relative to its value at the other receiving coil. The 
magnitude is often referred to as the attenuation, and the phase is often called the phase shift. 
Thus, the magnitude, or attenuation, and the phase shift, or phase, are typically derived from 
the ratio of the voltage at one receiver antenna relative to the voltage at another receiver 
antenna. 

[0007] Commercially deployed MWD/LWD resistivity measurement systems use multiple 
transmitters; consequently, attenuation and phase-based resistivity values can be derived 
independently using each transmitter or from averages of signals from two or more 
transmitters. See for example U.S. Patent 5,594,343. 

[0008] As demonstrated in U.S. Patents 4,968,940 and 4,899,1 12, a very common method 
practiced by those skilled in the art of MWD/LWD for determining the resistivity from the 
measured data is to transform the dielectric constant into a variable that depends on the 
resistivity and then to independently convert the phase shift and attenuation measurements to 
two separate resistivity values. A key assumption implicitly used in this practice is that each 
measurement senses the resistivity within the same volume that it senses the dielectric 
constant. This currently practiced method may provide significantly incorrect resistivity 
values, even in virtually homogeneous earth formations; and the errors may be even more 
severe in inhomogeneous formations. 

[0009] A MWD/LWD tool typically transmits a 2 MHz signal (although frequencies as low 
as 0.4 MHz are sometimes used). This frequency range is high enough to create difficulties 
in transforming the raw attenuation and phase measurements into accurate estimates of the 
resistivity and/or the dielectric constant. For example, the directly measured values are not 
linearly dependent on either the resistivity or the dielectric constant (this nonlinearity, known 
to those skilled in the art as "skin-effect," also limits the penetration of the fields into the 
earth formation). In addition, it is useful to separate the effects of the dielectric constant and 
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the resistivity on the attenuation and phase measurements given that both the resistivity and 
the dielectric constant typically vary spatially within the earth formation. If the effects of 
both of these variables on the measurements are not separated, the estimate of the resistivity 
can be corrupted by the dielectric constant, and the estimate of the dielectric constant can be 
corrupted by the resistivity. Essentially, the utility of separating the effects is to obtain 
estimates of one parameter that do not depend on (are independent of) the other parameter. A 
commonly used current practice relies on assuming a correlative relationship between the 
resistivity and dielectric constant (i.e., to transform the dielectric constant into a variable that 
depends on the resistivity) and then calculating resistivity values independently from the 
attenuation and phase shift measurements that are consistent the correlative relationship. 
Differences between the resistivity values derived from corresponding phase and attenuation 
measurements are then ascribed to spatial variations (inhomogeneities) in the resistivity over 
the sensitive volume of the phase shift and attenuation measurements. See for example U.S. 
Patents 4,899,1 12 and 4,968,940. An implicit and instrumental assumption in this method is 
that the attenuation measurement senses both the resistivity and dielectric constant within the 
same volume, and that the phase shift measurement senses both variables within the same 
volume (but not the same volume as the attenuation measurement). See for example U.S. 
Patents 4,899,112 and 4,968,940. These assumptions facilitate the independent 
determination of a resistivity value from a phase measurement and another resistivity value 
from an attenuation measurement. However, the implicit assumption mentioned above is not 
true; so, the results determined using such algorithms are questionable. 
[0010] Another technique for determining the resistivity and/or dielectric constant is to 
assume a model for the measurement apparatus in, for example, a homogeneous medium (no 
spatial variation in either the resistivity or dielectric constant) and then to determine values 
for the resistivity and dielectric constant that cause the model to agree with the measured 
phase shift and attenuation data. The resistivity and dielectric constant determined by the 
model are then correlated to the actual parameters of the earth formation. This method is 
thought to be valid only in a homogeneous medium because of the implicit assumption 
mentioned in the above paragraph. A recent publication by P.T. Wu, J.R. Lovell, B. Clark, 
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S.D. Bonner, and J.R. Tabanou entitled "Dielectric-Independent 2-MHz Propagation 
Resistivities" (SPE 56448, 1999) (hereafter referred to as "Wu") demonstrates that such 
assumptions are used by those skilled in the art. For example, Wu states: "One fundamental 
assumption in the computation of Rex is an uninvaded homogeneous formation. This is 
because the phase shift and attenuation investigate slightly different volumes." Abandoning 
the false assumptions applied in this practice results in estimates of one parameter (i.e., the 
resistivity or dielectric constant) that have no net sensitivity to the other parameter. This 
desirable and previously unknown property of the results is very useful because earth 
formations are commonly inhomogeneous. 

[0011] Wireline dielectric measurement tools commonly use electrical signals having 
frequencies in the range 20 MHz - 1.1 GHz. In this range, the skin-effect is even more 
severe, and it is even more useful to separate the effects of the dielectric constant and 
resistivity. Those skilled in the art of dielectric measurements have also falsely assumed that 
a measurement (either attenuation or phase) senses both the resistivity and dielectric constant 
within the same volume. The design of the measurement equipment and interpretation of the 
data both reflect this. See for example U.S. Patents 4,1 85,238 and 4,209,747. 
[0012] Wireline induction measurements are typically not attenuation and phase, but 
instead the real (R) and imaginary (X) parts of the voltage across a receiver antenna which 
consists of several induction coils in electrical series. For the purpose of this disclosure, the 
R-signal for a wireline induction measurement corresponds to the phase measurement of a 
MWD/LWD resistivity or wireline dielectric tool, and the X-signal for a wireline induction 
measurement corresponds to the attenuation measurement of a MWD/LWD resistivity or 
wireline dielectric device. Wireline induction tools typically operate using electrical signals 
at frequencies from 8-200 kHz (most commonly at approximately 20 kHz). This frequency 
range is too low for significant dielectric sensitivity in normally encountered cases; however, 
the skin-effect can corrupt the wireline induction measurements. As mentioned above, the 
skin-effect shows up as a non-linearity in the measurement as a function of the formation 
conductivity, and also as a dependence of the measurement sensitivity values on the 
formation conductivity. Estimates of the formation conductivity from wireline induction 
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devices are often derived from data processing algorithms which assume the tool response 
function is the same at all depths within the processing window. The techniques disclosed in 
the commonly-assigned U.S. Patent Application, Serial No. 09/877,383, entitled "Method of 
Determining Resistivity of an Earth Formation with Phase Resistivity Evaluation Based on a 
Phase Shift Measurement and Attenuation Resistivity Evaluation Based on an Attenuation 
Measurement and the Phase Shift Measurement," can be applied to wireline induction 
measurements for the purpose of deriving resistivity values without assuming the tool 
response function is the same at all depths within the processing window as is done in U.S. 
Patent 5,1 57,605. In order to make such an assumption, a background conductivity, a , that 
applies for the data within the processing window is commonly used. Practicing an 
embodiment as disclosed in the above referenced, commonly assigned U.S. Patent 
Application reduces the dependence of the results on the accuracy of the estimates for the 
background parameters because the background parameters are not required to be the same at 
all depths within the processing window. Another advantage is reduction of the need to 
perform steps to correct wireline induction data for the skin effect. Such parameterization as 
disclosed in the above referenced, commonly assigned U.S. Patent Application, however, has 
not fully taken into account complications such as invaded zones and frequency dispersion. 
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SUMMARY OF THE INVENTION 
[0013] Briefly, a method of estimating electrical parameters of an earth formation employs 
a simplified model of a measurement tool or device in transforming or normalizing data 
measured by the measurement tool. Electrical parameters of the earth formation, such as 
conductivity or dielectric constant, for example, may be estimated based on the normalized 
data. 

[0014] The foregoing has outlined rather broadly the features and technical advantages of 
the present invention in order that the detailed description of the invention that follows may 
be better understood. Additional features and advantages of the invention will be described 
hereinafter which form the subject of the claims of the invention. It should be appreciated by 
those skilled in the art that the conception and the specific embodiment disclosed may be 
readily utilized as a basis for modifying or designing other structures for carrying out the 
same purposes of the present invention. It should be also be realized by those skilled in the 
art that such equivalent constructions do not depart from the spirit and scope of the invention 
as set forth in the appended claims. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[0015] For a more complete understanding of the present invention, and the advantages 
thereof, reference is now made to the following descriptions taken in conjunction with the 
accompanying drawings, in which: 

[0016] FIGURE 1 is a plot of multiple laboratory measurements on rock samples 
representing the relationship between the conductivity and the dielectric constant in a variety 
of geological media; 

[0017] FIGURE 2 illustrates the derivation of a sensitivity function in relation to an 
exemplary one-transmitter, one-receiver MWD/LWD resistivity tool; 
[001 8] FIGURE 3 illustrates an exemplary one-transmitter, two-receiver MWD/LWD tool 
commonly referred to as an uncompensated measurement device; 

[0019] FIGURE 4 illustrates an exemplary two-transmitter, two-receiver MWD/LWD tool, 
commonly referred to as a compensated measurement device; 

[0020] FIGURES 5a, 5b, 5c and 5d are exemplary sensitivity function plots for Deep and 
Medium attenuation and phase shift measurements; 

[0021] FIGURES 6a, 6b, 6c and 6d are plots of the sensitivity functions for the Deep and 
Medium measurements of FIGURES 5a, 5b, 5c and 5d respectively transformed according to 
the techniques of a disclosed embodiment; 

[0022] FIGURE 7 is a portion of a table of background medium values and integral values 
employed in a disclosed embodiment; 

[0023] FIGURE 8 is a plot of attenuation and phase as a function of resistivity and 
dielectric constant; 

[0024] FIGURE 9 is a flowchart of a process that implements the techniques of a disclosed 
embodiment; 

[0025] FIGURE 10 is a table of exemplary data values for a medium spaced 2 MHz 
measurement described in conjunction with FIGURE 7; 

[0026] FIGURE 1 1 is a table of exemplary data values for a medium spaced 2 MHz 
measurement illustrating the difference between measured attenuation and phase-shift values 
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and corresponding point-dipole attenuation and phase-shift values for a fixed resistivity as a 
function of the dielectric constant; 

[0027] FIGURE 1 2 is a plot showing the difference between measured (mandrel) and point- 
dipole values as a function of both the resistivity and dielectric constant for the medium 
spaced 2 MHz measurement; and 

[0028] FIGURE 1 3 is a plot showing the performance of the extended approximations that 
facilitate use of rapidly-evaluated models. 
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DETAILED DESCRIPTION 
[0029] Some of the disclosed embodiments are relevant to both wireline induction and 
dielectric applications, as well as Measurement- While-Drilling and Logging-While-Drilling 
(MWD/LWD) applications. Turning now to the figures, Figure 1 is a plot of measurements 
of the conductivity and dielectric constant determined by laboratory measurements on a 
variety of rock samples from different geological environments. The points 121 through 129 
represent measured values of conductivity and dielectric constant (electrical parameters) for 
carbonate and sandstone earth formations. For instance, the point 126 represents a sample 
with a conductivity value of 0.01 (10" 2 ) Siemens per meter (S/m) and a relative dielectric 
constant of approximately 22. It should be noted that both the conductivity scale and the 
dielectric scale are logarithmic scales; so, the data would appear to be much more scattered if 
they were plotted on linear scales. 

[0030] The maximum boundary 1 1 1 indicates the maximum dielectric constant expected 
to be observed at each corresponding conductivity, hi a similar fashion, the minimum 
boundary 115 represents the minimum dielectric constant expected to be observed at each 
corresponding conductivity. Thepoints 122 through 128 represent measured values that fall 
somewhere in between the minimum boundary 115 and the maximum boundary 111. A 
median line 1 13 is a line drawn so that half the points, or points 121 through 124 are below 
the median line 113 and half thepoints, or points 126 through 129 are above the median line 
113. The point 125 falls right on top of the median line 113. 

[0031] An elemental measurement between a single transmitting coil 205 and a single 
receiving coil 207 is difficult to achieve in practice, but it is useful for describing the 
sensitivity of the measurement to variations of the conductivity and dielectric constant within 
a localized volume 225 of an earth formation 215. Figure 2 illustrates in more detail 
specifically what is meant by the term "sensitivity function," also referred to as a "response 
function" or "geometrical factors." Practitioners skilled in the art of wireline logging, 
Measurement- While-Drilling (MWD) and Logging- While-Drilling (LWD) are familiar with 
how to generalize the concept of a sensitivity function to apply to realistic measurements 
from devices using multiple transmitting and receiving antennas. Typically a MWD/LWD 
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resistivity measurement device transmits a signal using a transmitter coil and measures the 
phase and magnitude of the signal at one receiver antenna 307 relative to the values of the 
phase and the magnitude at another receiver antenna 309 within a borehole 301 (Fig. 3). 
These relative values are commonly referred to as the phase shift and attenuation. It should 
be understood that one way to represent a complex signal with multiple components is as a 
phasor signal. 

SENSITIVITY FUNCTIONS 

[0032] Figure 2 illustrates an exemplary single transmitter, single receiver MWD/LWD 
resistivity tool 220 for investigating an earth formation 215. A metal shaft, or "mandrel," 
203 is incorporated within the drill string (the drill string is not shown, but it is a series of 
pipes screwed together with a drill bit on the end), inserted into the borehole 201, and 
employed to take measurements of an electrical signal that originates at a transmitter 205 and 
is sensed at a receiver 207. The measurement tool is usually not removed from the well until 
the drill string is removed for the purpose of changing drill bits or because drilling is 
completed. Selected data from the tool are telemetered to the surface while drilling. All data 
are typically recorded in memory banks for retrieval after the tool is removed from the 
borehole 201 . Devices with a single transmitter and a single receiver are usually not used in 
practice, but they are helpful for developing concepts such as that of the sensitivity function. 
Schematic drawings of simple, practical apparatuses are shown in Figures 3 and 4. 
[0033] In a wireline operation, the measurement apparatus is connected to a cable (known 
as a wireline), lowered into the borehole 201, and data are acquired. This is done typically 
after the drilling operation is finished. Wireline induction tools measure the real (R) and 
imaginary (X) components of the receiver 207 signal. The R and X-signals correspond to the 
phase shift and attenuation measurements respectively. In order to correlate the sensitivity of 
the phase shift and attenuation measurements to variations in the conductivity and dielectric 
constant of the earth formation 215 at different positions within the earth formation, the 
conductivity and dielectric constant within a small volume P 225 are varied. For simplicity, 
the volume P 225 is a solid of revolution about the tool axis (such a volume is called a two- 
dimensional volume). The amount the phase and attenuation measurements change relative 
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to the amount the conductivity and dielectric constant changed within P 225 is essentially 
the sensitivity. The sensitivity function primarily depends on the location of the point P 225 
relative to the locations of the transmitter 205 and receiver 207, on the properties of the earth 
formation 2 1 5 , and on the excitation frequency. It also depends on other variables such as the 
diameter and composition of the mandrel 203, especially when P is near the surface of the 
mandrel 203. 

[0034] Although the analysis is carried out in two-dimensions, the important conclusions 
regarding the sensitive volumes of phase shift and attenuation measurements with respect to 
the conductivity and dielectric constant hold in three-dimensions. Consequently, the scope of 
this application is not limited to two-dimensional cases. This is discussed more in a 
subsequent section entitled, "ITERATIVE FORWARD MODELING AND DIPPING 
BEDS." 

[0035] The sensitivity function can be represented as a complex number having a real and 
an imaginary part. In the notation used below, S , denotes a complex sensitivity function, 
and its real part is S' , and its imaginary part is S" . Thus, S = S' + iS" , in which the 
imaginary number i = V^T - The quantities S' and S" are commonly referred to as 
geometrical factors or response functions. The volume P 225 is located a distance p in the 
radial direction from the tool's axis and a distance z in the axial direction from the receiver 
206. S' represents the sensitivity of attenuation to resistivity and the sensitivity of phase 
shift to dielectric constant. Likewise, S" represents the sensitivity of attenuation to dielectric 
constant and the sensitivity of phase shift to resistivity. The width of the volume P 225 is 
Ap 21 land the height of the volume P 225 is Az 213. The quantity S' , or the sensitivity of 
attenuation to resistivity, is calculated by determining the effect a change in the conductivity 
(reciprocal of resistivity) in volume P 225 from a prescribed background value has on the 
attenuation of a signal between the transmitter 205 and the receiver 207, assuming the 
background conductivity value is otherwise unperturbed within the entire earth formation 
215. In a similar fashion, S" , or the sensitivity ofthe phase to the resistivity, is calculated by 
determining the effect a change in the conductivity value in the volume P 225 from an 
assumed background conductivity value has on the phase of the signal between the 
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transmitter 205 and the receiver 207, assuming the background parameters are otherwise 
unperturbed within the earth formation 215. Alternatively, one could determine S' and 
S" by determining the effect a change in the dielectric constant within the volume P 225 has 
on the phase and attenuation, respectively. When the sensitivities are determined by 
considering a perturbation to the dielectric constant value within the volume P 225, it is 
apparent that the sensitivity of the attenuation to changes in the dielectric constant is the same 
as the sensitivity of the phase to the conductivity. It is also apparent that the sensitivity of the 
phase to the dielectric constant is the same as the sensitivity of the attenuation to the 
conductivity. By simultaneously considering the sensitivities of both the phase and 
attenuation measurement to the dielectric constant and to the conductivity, the Applicant 
shows a previously unknown relationship between the attenuation and phase shift 
measurements and the conductivity and dielectric constant values. By employing this 
previously unknown relationship, the Applicant provides techniques that produce better 
estimates of both the conductivity and the dielectric constant values than was previously 
available from those with skill in the art. The sensitivity functions S' and S" and their 
relation to the subject matter of the Applicant's disclosure is explained in more detail below 
in conjunction with Figures 5a-d and Figures 6a-d. 

[0036] In the above, sensitivities to the dielectric constant were referred to. Strictly 
speaking, the sensitivity to the radian frequency co times the dielectric constant should have 
been referred to. This distinction is trivial to those skilled in the art. 

[0037] In Figure 2, if the background conductivity (reciprocal of resistivity) of the earth 
formation 215 is a 0 and the background dielectric constant of the earth formation 215 is e 0 , 
then the ratio of the receiver 207 voltage to the transmitter 205 current in the background 
medium can be expressed as Z° r , where R stands for the receiver 207 and T stands for the 
transmitter 205. Hereafter, a numbered subscript or superscript such as the '0' is merely used 
to identify a specific incidence of the corresponding variable or function. If an exponent is 
used, the variable or function being raised to the power indicated by the exponent will be 
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surrounded by parentheses and the exponent will be placed outside the parentheses. For 
example (L,) 3 would represent the variable L x raised to the third power. 
[00381 When the background conductivity ct 0 and/or dielectric constant e 0 are replaced 
new values ct, and/or e x in the volume P 225, the ratio between the receiver 207 voltage to 
the transmitter 205 current is represented by Z\ T . Using the same nomenclature, a ratio 
between a voltage at a hypothetical receiver placed in the volume P 225 and the current at the 
transmitter 205 can be expressed as Z° PT . In addition, a ratio between the voltage at the 
receiver 207 and a current at a hypothetical transmitter in the volume P 225 can be expressed 
as Zip . Using the Born approximation, it can be shown that, 

^f- = 1 + S(T, R, P) AoApAz 
Z RT 

[0039] where the sensitivity function, defined as S(T,R,V), is 
S(T,R,P) = -f^ 

[0040] in which Act = ct, - ct 0 = (ct, - ct 0 ) + 1 co(e x - e 0 ) , and the radian frequency of the 
transmitter current is « = 2^f. A measurement of this type, in which there is just one 
transmitter 205 and one receiver 207, is defined as an "elemental" measurement. It should be 
noted that the above result is also valid if the background medium parameters vary spatially 
within the earth formation 215. In the above equations, both the sensitivity function 
S(T,R,P) and the perturbation Act are complex-valued. Some disclosed embodiments 
consistently treat the measurements, their sensitivities, and the parameters to be estimated as 
complex-valued functions. This is not done in the prior art. 

[0041 ] The above sensitivity function of the form S(T, R, P) is referred to as a 2-D (or 
two-dimensional) sensitivity function because the volume ApAz surrounding the point P 
225, is a solid of revolution about the axis of the tool 201 . Because the Born approximation 
was used, the sensitivity function S depends only on the properties of the background medium 



Docket No. PAT009CON PATENT 

16 

because it is assumed that the same field is incident on the point P(p,z) even though the 
background parameters have been replaced by <r, . 

[0042] Figure 3 illustrates an exemplary one-transmitter, two-receiver MWD/LWD 
resistivity measurement apparatus 320 for investigating an earth formation 315. Due to its 
configuration, the tool 320 is defined as an "uncompensated" device and collects 
uncompensated measurements from the earth formation 315. For the sake of simplicity, a 
borehole is not shown. This measurement tool 320 includes a transmitter 305 and two 
receivers 307 and 309, each of which is incorporated into a metal mandrel 303. Typically, 
the measurement made by such a device is the ratio of the voltages at receivers 307 and 309. 
In this example, using the notation described above in conjunction with Figure 2, the 
sensitivity function S(T,R,R\P) for the uncompensated device can be shown to be the 
difference between the elemental sensitivity functions S(T,R,P) and S(T,R',P), where T 
represents the transmitter 305, R represents the receiver 307, R' represents the receiver 309, 
and P represents a volume (not shown) similar to the volume P 225 of Figure 2. 
[0043] For wireline induction measurements, the voltage at the receiver R is subtracted 
from the voltage at the receiver R' , and the position and number of turns of wire for R are 
commonly chosen so that the difference in the voltages at the two receiver antennas is zero 
when the tool is in a nonconductive medium. For MWD/LWD resistivity and wireline 
dielectric constant measurements, the voltage at the receiver R , or V R , and the voltage at the 
receiver R' , or V R , , are examined as the ratio V R jV R , . In either case, it can be shown that 

S(T, R, R', P) = S{T,R,P) - S(T, R', P). 
[0044] The sensitivity for an uncompensated measurement is the difference between the 
sensitivities of two elemental measurements such as S(T,R,P) and S(T,R',P) calculated as 
described above in conjunction with Figure 2. 

[0045] Figure 4 illustrates an exemplary two-transmitter, two-receiver MWD/LWD 
resistivity tool 420. Due to its configuration (transmitters being disposed symmetrically), the 
tool 420 is defined as a "compensated" tool and collects compensated measurements from an 
earth formation 415. The tool 420 includes two transmitters 405 and 41 1 and two receivers 
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407 and 409, each of which is incorporated into a metal mandrel or collar 403. Each 
compensated measurement is the geometric mean of two corresponding uncompensated 
measurements. In other words, during a particular timeframe, the tool 420 performs two 
uncompensated measurements, one employing transmitter 405 and the receivers 407 and 409 
and the other employing the transmitter 411 and the receivers 409 and 407. These two 
uncompensated measurements are similar to the uncompensated measurement described 
above in conjunction with Figure 3. The sensitivity function S of the tool 420 is then defined 
as the arithmetic average of the sensitivity functions for each of the uncompensated 
measurements. Another way to describe this relationship is with the following formula: 

S(T,R, R', T' t P) = ~ [S(T, R,R',P) + S(T', R',R, P)] 
[0046] where T represents transmitter 405 , V represents transmitter 41 1 , R represents 
receiver 407, R' represents receiver 409 and P represents a small volume of the earth 
formation similar to 225 (Fig. 2). 

[0047] The techniques of the disclosed embodiments are explained in terms of a 
compensated tool such as the tool 420 and compensated measurements such as those 
described in conjunction with Figure 4. However, it should be understood that the techniques 
also apply to uncompensated tools such as the tool 320 and uncompensated measurements 
described above in conjunction with Figure 3 and elemental tools such as the tool 220 and 
elemental measurements such as those described above in conjunction with Figure 2. In 
addition, the techniques are applicable for use in a wireline system, a system that may not 
incorporate its transmitters and receivers into a metal mandrel, but may rather affix a 
transmitter and a receiver to a tool made of a non-conducting material such as fiberglass. The 
wireline induction frequency is typically too low for dielectric effects to be significant. Also 
typical for wireline induction systems is to select the position and number of turns of groups 
of receiver antennas so that there is a null signal in a nonconductive medium. When this is 
done, Z° RT = 0 if <r 0 = 0 . As a result, it is necessary to multiply the sensitivity and other 
quantities by Z° RT to use the formulation given here in such cases. 
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[0048] The quantity Z\ T /Z° RT can be expressed as a complex number which has both a 
magnitude and a phase (or alternatively real and imaginary parts). To a good approximation, 
the raw attenuation value (which corresponds to the magnitude) is: 



-g 2 - = 1 + Re[S(r, R, P)Aa]ApAz = 1 + [S'Acx - S"caA£]ApAz 

z RT 

[0049] where the function Re[ ] denotes the real part of its argument. Also, to a good 
approximation, the raw phase shift value is: 

phas J^M. j » Imjj^J = \m[S(T,R, P)Aa]ApAz = [S"Acr + S'coAs ]ApAz 

[0050] in which Im[] denotes the imaginary part of its argument, S(T,R,P) = S' + iS" , 
Act = <7, - ct 0 , and As = e l -e 0 . For the attenuation measurement, S' is the sensitivity to 
the resistivity and S" is the sensitivity to the dielectric constant. For the phase shift 
measurement, S' is the sensitivity to the dielectric constant and S" is the sensitivity to the 
resistivity. This is apparent because S' is the coefficient of Act in the equation for 
attenuation, and it is also the coefficient of coAs in the equation for the phase shift. 
Similarly, S" is the coefficient of Act in the equation for the phase shift, and it is also the 
coefficient for - coAe in the equation for attenuation. This implies that the attenuation 
measurement senses the resistivity in the same volume as the phase shift measurement senses 
the dielectric constant and that the phase shift measurement senses the resistivity in the same 
volume as the attenuation measurement senses the dielectric constant. In the above, we have 
referred to sensitivities to the dielectric constant. Strictly speaking, the sensitivity to the 
radian frequency co times the dielectric constant Ae should have been referred to. This 
distinction is trivial to those skilled in the art. 

[0051] The above conclusion regarding the volumes in which phase and attenuation 
measurements sense the resistivity and dielectric constant from Applicant's derived equations 
also follows from a well known result from complex variable theory known in that art as the 
Cauchy-Reimann equations. These equations provide the relationship between the 
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derivatives of the real and imaginary parts of an analytic complex function with respect to the 
real and imaginary parts of the function's argument. 

[00521 Figures 5a, 5b, 5c and 5d can best be described and understood together. In all 
cases, the mandrel diameter is 6.75 inches, the transmitter frequency is 2 MHz, and the 
background medium is characterized by a conductivity of <r 0 = 0.01 S/m and a relative 
dielectric constant of =10. The data in Figures 5a and 5c labeled "Medium 
Measurement" are for a compensated type of design shown in Figure 4. The exemplary 
distances between transmitter 405 and receivers 407 and 409 are 20 and 30 inches, 
respectively. Since the tool is symmetric, the distances between transmitter 41 1 andrecervers 
409 and 407 are 20 and 30 inches, respectively. The data in Figures 5b and 5d labeled "Deep 
Measurement" are also for a compensated tool as shown in Figure 4, but with exemplary 
transmitter-receiver spacings of 50 and 60 inches. Each plot shows the sensitivity of a given 
measurement as a function of position within the formation. The term sensitive volume 
refers to the shape of each plot as well as its value at any point in the formation. The axes 
labeled"Axial Distance" refer to the coordmate along the axis of the tool with zero beingthe 
geometric mid-point of the antenna array (halfway between receivers 407 and 409) to a given 
point in the formation. The axes labeled "Radial Distance" refer to the radial distance from 
theaxisofthetooltoagivenpointintheformation. The value on the vertical axis is actually 
the sensitivity value for the indicated measurement. Thus, Figure 5ais a plot of a sensitivity 
function that illustrates the sensitivity ofthe "Medium" phase shift measurement in relation 
to changes in the resistivity as a function of the location of the point P 225 in the earth 
formation 21 5 (Fig. 2). If the measurement of phase shift changes significantly in response to 
changing the resistivity from its background value, then phase shift is considered relatively 
sensitive to the resistivity at the point P 225. If the measurement of phase shift does not 
change significantly in response to changing the resistivity, then the phase sWft is con^ 
relatively insensitive at the point P 225 . Based upon the relationship disclosed herein, Figure 
5a also illustrates the sensitivity ofthe "Medium" attenuation measurement in relation to 
changes in dielectric constant values. Note that the dimensions ofthe sensitivity on the 
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vertical axes is ohms per meter (Q/m) and distances on the horizontal axes are listed in 
inches. In a similar fashion, Figure 5b is a plot of the sensitivity of the attenuation 
measurement to the resistivity. Based on the relationship disclosed herein, Figure 5b is also 
the sensitivity of a phase shift measurement to a change in the dielectric constant. Figures 5b 
and 5d have the same descriptions as Figures 5a and 5c, respectively, but Figures 5b and 5d 
are for the "Deep Measurement" with the antenna spacings described above. 
[0053] Note that the shape of Figure 5a is very dissimilar to the shape of Figure 5c. This 
means that the underlying measurements are sensitive to the variables in different volumes. 
For example, the Medium phase shift measurement has a sensitive volume characterized by 
Figure 5 a for the resistivity, but this measurement has the sensitive volume shown in Figure 
5c for the dielectric constant. Figures 5a-5d illustrate that for a particular measurement the 
surface S" is more localized than the surface S' such that the sensitive volume associated 
with the surface S' substantially encloses the sensitive volume associated with S" for the 
corresponding measurement. As discussed below, it is possible to transform an attenuation 
and a phase shift measurement to a complex number which has the following desirable 
properties: 1) its real part is sensitive to the resistivity in the same volume that the imaginary 
part is sensitive to the dielectric constant; 2) the real part has no net sensitivity to the 
dielectric constant; and, 3) the imaginary part has no net sensitivity to the resistivity. In 
addition, the transformation is generalized to accommodate multiple measurements acquired 
at multiple depths. The generalized method can be used to produce independent estimates of 
the resistivity and dielectric constant within a plurality of volumes within the earth formation. 
TRANSFORMED SENSITIVITY FUNCTIONS AND TRANSFORMATION OF THE 
MEASUREMENTS 

[0054] For simplicity, the phase shift and attenuation will not be used. Hereafter, the real 
and imaginary parts of measurement will be referred to instead. Thus, 




w'=(10) d 
w"=(10)' 



x cos(6>) 
' x sin(#) 
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[0055] where w' is the real part of w, w" is the imaginary part of w, i is the square root 
of the integer -1, dB is the attenuation in decibels, and 0 is the phase shift in radians. 
[0056[ The equations that follow can be related to the sensitivity functions described 
above in conjunction with Figure 2 by defining variables w, = Z\ T and w 0 = Z° RT . The 
variable w, denotes an actual tool measurement in the earth formation 215. The variable w 0 
denotes the expected value for the tool measurement in the background earth formation 215. 
For realistic measurement devices such as those described in Figures 3 and 4, the values for 
w, and w 0 would be the voltage ratios defined in the detailed description of Figures 3 and 4. 
In one embodiment, the parameters for the background medium are determined and then used 
to calculate value of w 0 using a mathematical model to evaluate the tool response in the 
background medium. One of many alternative methods to determine the background medium 
parameters is to estimate w 0 directly from the measurements, and then to determine the 
background parameters by correlating w 0 to a model of the tool in the formation which has 
the background parameters as inputs. 

[0057] As explained in conjunction with Figure 2, the sensitivity function S relates the 
change in the measurement to a change in the medium parameters such as resistivity and 
dielectric constant within a small volume 225 of the earth formation 215 at a prescribed 
location in the earth formation 225, or background medium. A change in measurements due 
to small variations in the medium parameters at a range of locations can be calculated by 
integrating the responses from each such volume in the earth formation 215. Thus, if Act is 
defined for a large number of points p, z , then 

Z RT =Z° RT (l + I[SAa]) 
[0058] in which / is a spatial integral function further defined as 

I[F]= ]dz]dpF(p,z) 

-» 0 

[0059] where F is a complex function. 
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[0060] Although the perturbation from the background medium, Act is a function of 
position, parameters of a hypothetical, equivalent homogeneous perturbation (meaning that 
no spatial variations are assumed in the difference between the resistivity and dielectric 
constant and values for both of these parameters in the background medium) can be 
determined by assuming the perturbation is not a function of position and then solving for it. 
Thus, 



[0061 ] where Act represents the parameters of the equivalent homogeneous perturbation. 
From the previous equations, it is clear that 



[0063] where Act is the transformed measurement (it is understood that Act is also the 
equivalent homogeneous perturbation and that the terms transformed measurement and 



function for the transformed measurement, and S will be referred to as the transformed 
sensitivity function. In the above, w, is the actual measurement, and w 0 is the value 
assumed by the measurement in the background medium. An analysis of the transformed 
sensitivity function S, shows that the transformed measurements have the following 
properties: 1) the real part of Act is sensitive to the resistivity in the same volume that its 
imaginary part is sensitive to the dielectric constant; 2) the real part of Act has no net 
sensitivity to the dielectric constant; and, 3) the imaginary part of Act has no net sensitivity 
to the resistivity. Details of this analysis will be given in the next few paragraphs. 
[0064] The techniques of the disclosed embodiment can be further refined by introducing 
a calibration factor c (which is generally a complex number that may depend on the 
temperature of the measurement apparatus and other environmental variables) to adjust for 



A6I[S] = I[SAcr] 




[0062] 



and 




equivalent homogeneous perturbation will be used synonymously), S is the sensitivity 
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anomalies in the physical measurement apparatus. In addition, the term, w bh can be 
introduced to adjust for effects caused by the borehole 201 on the measurement. With these 
modifications, the transformation equation becomes 



Act = = USA*] = J-f^^ -ll. 

I[S] I[S]{ w 0 ) 



[0065] The sensitivity function for the transformed measurement is determined by 
applying the transformation to the original sensitivity function, S . Thus, 
c_^^- S S'I[S'] + S"I[S"} S"I[S']-S'I[S"] 

s = s+lS= i[S] = \ J[S f \i[sf ■ 

[0066] Note that I[S] = I[S'] = 1 because I[S"] = 0. The parameters for the equivalent 
homogeneous perturbation are 

Act' = cr, - cr 0 = I[S'Aa] - I[S"o)As] 

Act" = co{s, -s 0 ) = I[S'o)Ae] + I[S"Act]. 
[0067] The estimate for the conductivity perturbation, Act' suppresses sensitivity (is 
relatively insensitive) to the dielectric constant perturbation, and the estimate of the dielectric 
constant perturbation, Ad"/co suppresses sensitivity to the conductivity perturbation. This is 
apparent because the coefficient of the suppressed variable is S" . In fact, the estimate for the 
conductivity perturbation Act ' is independent of the dielectric constant perturbation provided 
that deviations in the dielectric constant from its background are such that J[S"(oA£ ] = 0 . 
Since I[S"] = 0 , this is apparently the case if coAe is independent of position. Likewise, the 
estimate for the dielectric constant perturbation given by Aa" / co is independent of the 
conductivity perturbation provided that deviations in the conductivity from its background 
value are such that I[S"Aa] = 0. Since I[S"] = 0, this is apparently the case if Act is 
independent of position. 

[0068] Turning now to Figures 6a and 6b, illustrated are plots of the sensitivity functions 
S' and S" derived from S' and S" for the medium transmitter-receiver spacing 
measurement shown in Figures 5a and 5c using the transformation 
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I[S] 



[0069] The data in Figures 6c and 6d were derived from the data in Figures 5b and 5d for 
the Deep T-R spacing measurement. As shown in Figures 6a, 6b, 6c and 6d, using the 
transformed measurements to determine the electrical parameters of the earth formation is a 
substantial improvement over the prior art. The estimates of the medium parameters are 
more accurate and less susceptible to errors in the estimate of the background medium 
because the calculation of the resistivity is relatively unaffected by the dielectric constant and 
the calculation of the dielectric constant is relatively unaffected by the resistivity. In addition 
to integrating to 0, the peak values for S" in Figures 6b and 6d are significantly less than the 
respective peak values for S' in Figures 6a and 6c. Both of these properties are very 
desirable because S" is the sensitivity function for the variable that is suppressed. 
REALIZATION OF THE TRANSFORMATION 

[0070] In order to realize the transformation, it is desirable to have values of I[S] readily 

accessible over the range of background medium parameters that will be encountered. One 

way to achieve this is to compute the values for J[S] and then store them in a lookup table 

for use later. Of course, it is not necessary to store these data in such a lookup table if it is 

practical to quickly calculate the values for I[S] on command when they are needed. In 

general, the values for I[S] can be computed by directly; however, it can be shown that 

1 dw I 
I[S] = — — 

W 0 d(T\5=a 0 

[0071] where w 0 is the expected value for the measurement in the background medium, 
and the indicated derivative is calculated using the following definition: 

dw \ ,• w(g 0 +Aa) -w(go) 
= lim — — — . 

da\s=5 Q a^o (cr 0 + Aa)-o- 0 

[0072] In the above formula, a 0 may vary from point to point in the formation 2 1 5 (the 
background medium may be inhomogeneous), but the perturbation Ac? is constant at all 
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points in the formation 215. As an example of evaluating I[S] using the above formula, 
consider the idealized case of a homogeneous medium with a small transmitter coil and two 
receiver coils spaced a distance L x and L 2 from the transmitter. Then, 



w n = 



LA cxp(ik 0 L 2 )(l-ik 0 L 2 ) 
L 2 ) expO^Xl-^oA) 
1 dw I ico M ( {L 2 f (A) 2 



I\S~\ — 

w 0 da\a=a 0 2 {l-ik Q L 2 \-ik 0 L x 



[0073] The wave number in the background medium is k 0 = ^iona Q , the function 
exp()is the complex exponential function where exp(l) » 2.71828 , and the symbol 
denotes the magnetic permeability of the earth formation. The above formula for I[S] 
applies to both uncompensated (Figure 3) and to compensated (Figure 4) measurements 
because the background medium has reflection symmetry about the center of the antenna 
array in Figure 4. 

[0074] For the purpose of this example, the above formula is used to compute the values 
for I[S] = I[S'] + HIS"] . Figure 7 illustrates an exemplary table 701 employed in a Create 
Lookup Table step 903 (Fig. 9) of the technique of the disclosed embodiment. Step 903 
generates a table such as table 701 including values for the integral of the sensitivity function 
over the range of variables of interest. The first two columns of the table 701 represent the 
conductivity cr 0 and the dielectric constant s 0 of the background medium. The third and 
fourth columns of the table 701 represent calculated values for the functions I[S'] and 
I[S n ] for a Deep measurement, in which the spacing between the transmitter 305 receivers 
307 and 309 is 50 and 60 inches, respectively. The fifth and sixth columns of the table 701 
represent calculated values for the functions I[S'] and I[S"] for a Medium measurement, in 
which the spacing between the transmitter 305 receivers 307 and 309 is 20 and 30 inches, 
respectively. It is understood that both the frequency of the transmitter(s) and the spacing 
between the transmitter(s) and receiver(s) can be varied. Based upon this disclosure, it is 
readily apparent to those skilled in the art that algorithms such as the one described above can 
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be applied to alternative measurement configurations. If more complicated background media 
are used, for example including the mandrel with finite-diameter antennas, it may be more 
practical to form a large lookup table such as table 701 bu, with many more values. Instead 
of calculating I[S] every time a value is needed, data would be interpolated from to table. 
Nonetheless, table 701 clearly illustrates the nature of such a lookup table. Such a table 
would contain the valuesofthe functions /[S'l and IIS'] for the entire range of values of 
the conductivity cr 0 and the dielectric constant e, likely to be encountered in typical earth 
formations. For example, f[S'] and US'] could be calculated for values of , 0 between 1 
and 1000 and for values of a 0 between 0.0001 and 10.0. Whether calculating values forthe 
entire looknp table 701 or eomputingthe I[S'] and /IS"] on command as needed, the data is 
used as explained below. 

[0075] Figure 8 illustrates a chart 801 used to tmplement a Determine Background 
Medium Parameters step 905 (Fig. 9) of the techniques of the drsclosed embodiment. The 
chart 801 represents a plot of the attenuation and phase shift as a function of resistivny and 
dielectric constant in a homogeneous medium. Similar plots can be derived for more 
complicated media. However, the homogeneous background media are routinelyused due to 
their simplicity. Well known numerical methods such as inverse interpolation can be used to 
calculate an initial estimate of background parameters based upon the chart 801. In one 
embodiment, the measured attenuation and phase shift values are averaged over a few feet of 
depth within the borehole 201 . These average values are used to determine the background 
resistivity and dielectric constant based upon the chart 801. 1. should be understood that 
background medium parameters can be estimated in a variety of ways using one or more 
attenuation and phase measurements. 

[0076] Figure 9 is a flowchart of an embodiment of the disclosed transformation 
techniques mat can be implemented in a software program which is executed by aprocessor 
of a computing system such as a computer a. the surface or a "downhole" microprocessor. 
Starting in a Begin Analysis step 901, control proceeds immediately to the Create lookup 
Table step 903 described above in conjunction with Figure 7. In an alternative embodrment, 
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step 903 can be bypassed and the function of the lookup table replaced by curve matching, or 
"forward modeling." Control then proceeds to an Acquire Measured Data Step 904. Next, 
control proceeds to a Determine Background Values Step 905, in which the background 
values for the background medium are determined. Step 905 corresponds to the chart 801 
(Fig. 8). 

[0077] Control then proceeds to a Determine Integral Value step 907. The Determine 
Integral Value step 907 of the disclosed embodiment determines an appropriate value for 
I[S] using the lookup table generated in the step 903 described above or by directly 
calculating the I[S] value as described in conjunction with Figure 7. Compute Parameter 
Estimate, step 909, computes an estimate for the conductivity and dielectric constant as 
described above using the following equation: 

[0078] where the borehole effect and a calibration factor are taken into account using the 
factors w bh and c , respectively. The conductivity value plotted on the log (this is the value 
correlated to the conductivity of the actual earth formation) is Re(Aa+a 0 ) where the 
backgroundmediumischaracterizedby a 0 . The estimate for the dielectric constant can also 
be plotted on the log (this value is correlated to the dielectric constant of the earth formation), 
and this value is Im(Acr + <r 0 )/» • Lastly, in the Final Depth step 911, it is determined 
whether the tool 201 is at the final depth within the earth formation 215 that will be 
considered in the current logging pass. If the answer is "Yes," then control proceeds to a step 
921 where is processing is complete. If the answer in step 91 1 is "No," control proceeds to a 
Increment Depth step 913 where the tool 220 is moved to its next position in the borehole 
201 which penetrates the earth formation 215. After incrementing the depth of the tool 220, 
control proceeds to step 904 where the process of steps 904, 905, 907, 909 and 911 are 
repeated. It should be understood by those skilled in the art that embodiments descnbed 
herein in the form a computing system or as a programmed electrical circuit can be realized. 
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[0079] Improved estimates for the conductivity and/or dielectric constant can be 
determined by simultaneously considering multiple measurements at multiple depths. This 
procedure is described in more detail below under the heading "Multiple Sensors At Multiple 
Depths." 

MULTIPLE SENSORS AT MULTIPLE DEPTHS 

[0080] In the embodiments described above, the simplifying assumption that Aa is not 
position dependent facilitates determining a value for Act associated with each measurement 
by considering only that measurement at a single depth within the well (at least given a 
background value a 0 ). It is possible to eliminate the assumption that Acr is independent of 
position by considering data at multiple depths, and in general, to also consider multiple 
measurements at each depth. An embodiment of such a technique for jointly transforming 
data from multiple MWD/LWD sensors at multiple depths is given below. Such an 
embodiment can also be used for processing data from a wireline dielectric tool or a wireline 
induction tool. Alternate embodiments can be developed based on the teachings of this 
disclosure by those skilled in the art. 

[0081] In the disclosed example, the background medium is not assumed to be the same 
at all depths within the processing window. In cases where it is possible to assume the 
background medium is the same at all depths within the processing window, the system of 
equations to be solved is in the form of a convolution. The solution to such systems of 
equations can be expressed as a weighted sum of the measurements, and the weights can be 
determined using standard numerical methods. Such means are known to those skilled in the 
art, and are referred to as "deconvolution" techniques. It will be readily understood by those 
with skill in the art that deconvolution techniques can be practiced in conjunction with the 
disclosed embodiments without departing from the spirit of the invention, but that the 
attendant assumptions are not necessary to practice the disclosed embodiments in general. 
[0082] Devices operating at multiple frequencies are considered below, but 
multifrequency operation is not necessary to practice the disclosed embodiments. Due to 
frequency dispersion (i.e., frequency dependence of the dielectric constant and/or the 
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conductivity value), it is not necessarily preferable to operate using multiple frequencies. 
Given the disclosed embodiments, it is actually possible to determine the dielectric constant 
and resistivity from single-frequency data. In fact, the disclosed embodiments can be used to 
determine and quantify dispersion by separately processing data sets acquired at different 
frequencies. In the below discussion, it is understood that subsets of data from a given 
measurement apparatus or even from several apparatuses can be processed independently to 
determine parameters of interest. The below disclosed embodiment is based on using all the 
data available strictly for purpose of simplifying the discussion. 

[0083] Suppose multiple transmitter-receiver spacings are used and that each transmitter 
is excited using one or more frequencies. Further, suppose data are collected at multiple 
depths in the earth formation 215. Let TV denote the number of independent measurements 
performed at each of several depths, where a measurement is defined as the data acquired at a 
particular frequency from a particular set of transmitters and receivers as shown in Figures 3 
or 4. Then, at each depth z k , a vector of all the measurements can be defined as 



A, K ) k2 V w o ) u 



[0084] and the perturbation of the medium parameters from the background medium 
values associated with these measurements is 

A^ = Aa(p ) z)[l,l,...,l] r 
[0085] in which the superscript T denotes a matrix transpose, v k is a vector each element 
of which is a measurement, and is a vector each element of which is a perturbation from 
the background medium associated with a corresponding element of v k at the point P 225. In 
the above, the dependence of the perturbation, A^(p, z) on the position of the point P 225 is 
explicitly denoted by the variables p and z . In general, the conductivity and dielectric 
constant of both the background medium and the perturbed medium depend on p and z ; 
consequently, no subscript k needs to be associated with Acr(p, z) , and all elements of the 
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vector A* are equal. As descnbed above, borehole corrections and a calibration can be 
applied to each measurement, but here they are omitted for simplicity. 
[0086] The vectors v k and Ac? are related as follows: 
v k = 7[I AcF] 

■0087, inwhich I isadiagonal,natnxwi,heachdiago„aU.e m e„tbein g «hese„si.ivity 
function centered on .he depth „ , for the corresponding dement of v.. and the m.egral 
operator I is defined by: 

/[F]= ]dz ]dpF(p,z). 

-oo 0 

[0088] Using the notation 

I mn [F]= \dz \dpF(p,z) 
l0 089, to denote integrals of a function over the indicated limns of integration, it is 
apparent that 



[0090, if P„ =0, p, = + », x*. — 

an integra, equation from which an estimate of Aff<p.z) can be calculated. W.th the 
definitions 

,0091, p:=(P. + P.-,)/2 and and making the approbation 

Aff(p,z) - A*0»:.O within the volumes associated with each value for m and n , it 
follows that 

[0092, where W ' < W to ensure this system of equations is not underdetermined. The 
unknown values Ao^O can then he determined hy soiving me above set of linear 
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equations. It is apparent that the embodiment described in the section entitled 
"REALIZATION OF THE TRANSFORMATION" is a special case of the above for which 
M = 0,N = N' = \. 

[0093] Although the approximation Aa(p, z) = Aa(p* , z* ) (which merely states that 
Aff(p,z) is a piecewise constant function of p,z) is used in the immediately above 
embodiment, such an approximation is not necessary. More generally, it is possible to 
expand Aa(p, z) using a set of basis functions, and to then solve the ensuing set of equations 
for the coefficients of the expansion. Specifically, suppose 

Aa(p,z) = J J a m „<P m „(P> z ) 

[0094] then, 

[0095] where a mn = a mn [\X...,\f . Some desirable properties for the basis functions <f> mn 
are: 1) the integrals l(E<j> mn ] in the above equation all exist; and, 2) the system of equations 
for the coefficients a mn is not singular. It is helpful to select the basis functions so that a 
minimal number of terms is needed to form an accurate approximation to AS{p,z) . 
[0096] The above embodiment is a special case for which the basis functions are unit step 
functions. In fact, employing the expansion 

Aa(p,z) = 2 ZA^p^^XM^-zJ-MCz-z^^twCp-pJ-^P-P-.)] 

m=-M n=l 

[0097] where «(•) denotes the unit step function leads directly to the same system of 
equations 

m=-M «=! 

[0098] given in the above embodiment. Specific values for M , N', z m , and p„ needed to 
realize this embodiment of the invention depend on the excitation frequency(ies), on the 
transmitter-receiver spacings that are under consideration, and generally on the background 
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conductivity and dielectric constant. Different values for z m and p n are generally used for 
different depth intervals within the same well because the background medium parameters 
vary as a function of depth in the well. 

[0099] Solving the immediately above system of equations results in estimates of the 
average conductivity and dielectric constant within the volume ofthe earth formation 215 
corresponding to each integral I mn [S ] . In an embodiment, the Least Mean Square method is 
used to determine values for ^{p n ,z m ) by solving the above system of equations. Many 
texts on linear algebra list other techniques that may also be used. 

[00100] Unlike other procedures previously used for processing MWD/LWD data, the 
techniques of a disclosed embodiment account for dielectric effects and provide for radial 
inhomogeneities in addition to bedding interfaces by consistently treating the signal as a 
complex-valued function ofthe conductivity and the dielectric constant. This procedure 
produces estimates of one variable (i.e., the conductivity) are not corrupted by effects ofthe 
other (i.e., the dielectric constant). 

[00101] A series of steps, similar to those of Figure 9, can be employed in order to 
implement the embodiment for Multiple Sensors at Multiple Depths. Since the lookup table 
for 7 m „[f ] needed to realize such an embodiment could be extremely large, these values are 
evaluated as needed in this embodiment. This can be done in a manner analogous to the 
means described in the above section "REALIZATION OF THE TRANSFORMATION" 
using the following formulae: 

1 dw I 
I nm [S] =— -==-\ 



dw I _ lim w(a 0 + Aa mn )-w(o- 0 ) 
aa mn |c=5 0 ~ A5 ^° (5 0 +A5 nm )-a 0 
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[00102] where S mn = a mn + icos mn represents the conductivity and dielectric constant of 
the region of space over which the integral I mn [S] is evaluated. In words, I mn [S] can be 
calculated by evaluating the derivative of the measurement with respect to the medium 
parameters within the volume covered by the integration. Alternatively, one could evaluate 

I mn [S] by directly carrying out the integration as needed. This eliminates the need to store 

the values in a lookup table. 

[001031 While the above exemplary systems are described in the context of an 
MWD/LWD system, it shall be understood that a system according to the described 
techniques can be implemented in a variety of other logging systems such as wireline 
induction or wireline dielectric measurement systems. Further in accordance with the 
disclosed techniques, it should be understood that phase shift and attenuation can be 
combined in a variety of ways to produce a component sensitive to resistivity and relatively 
insensitive to dielectric constant and a component sensitive to dxelectric constant and 
relatively insensitive to resistivity. In the instance of MWD/LWD resistivity measurement 
systems, resistivity is the variable of primary interest; as a result, phase shift and attenuation 
measurements canbe combined to produce a component sensitive to resistivity and relatively 
insensitive dielectric constant. 

SINGLE MEASUREMENTS AT A SINGLE DEPTH 

[001 04] One useful embodiment is to correlate (or alternatively equate) a single measured 
value w, to a model that predicts the value of the measurement as a function of the 
conductivity and dielectric constant within a prescribed region of the earth formation. The 
value for the dielectric constant and conductivity that provides an acceptable correlation (or 
alternatively solves the equation) is then used as the final result (i.e., correlated to the 
parameters of the earth formation). This procedure can be performed mathematically, or 
graphically. Plotting a point on a chart such as Figure 8 and then determining which 
dielectric value and conductivity correspond to it is an example of performing the procedure 
graphically. It can be concluded from the preceding sections, that S is the sensitivity of such 
an estimate of the dielectric constant and conductivity to perturbations in either variable. 
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Thus such a procedure results in an estimate for the conductivity that has no net sensitivity to 
the changes in the dielectric constant and an estimate for the dielectric constant that has no 
net sensitivity to changes in the conductivity within the volume in question. This is a very 
desirable property for the results to have. The utility of employing a single measurement at a 
single depth derives from the fact that data processing algorithms using minimal data as 
inputs tend to provide results quickly and reliably. This procedure is a novel means of 
determining one parameter (either the conductivity or the dielectric constant) with no net 
sensitivity to the other parameter. Under the old assumptions, this procedure would appear to 
not be useful for determining independent parameter estimates. 
ITERATIVE FORWARD MODELING AND DIPPING BEDS 

[00105] The analysis presented above has been carried out assuming a 2-dimensional 
geometry where the volume P 225 in Figure 2 is a solid of revolution about the axis of the 
tool. In MWD/LWD and wireline operations, there are many applications where such a 2- 
dimensional geometry is inappropriate. For example, the axis of the tool often intersects 
boundaries between different geological strata at an oblique angle. Practitioners refer to the 
angle between the tool axis and a vector normal to the strata as the relative dip angle. When 
the relative dip angle is not zero, the problem is no longer 2-dimensional. However, the 
conclusion that: 1) the attenuation measurement is sensitive to the conductivity in the same 
volume as the phase measurement is sensitive to the dielectric constant; and, 2) an attenuation 
measurement is sensitive to the dielectric constant in the same volume that the phase 
measurement is sensitive to the conductivity remains true in the more complicated geometry. 
Mathematically,'this conclusion follows from the Cauchy-Reimann equations which still 
apply in the more complicated geometry (see the section entitled "SENSITIVITY 
FUNCTIONS"). The physical basis for this conclusion is that the conduction currents are in 
quadrature (90 degrees out of phase) with the displacement currents. At any point in the 
formation, the conduction currents are proportional to the conductivity and the displacement 
currents are proportional to the dielectric constant. 

[001061 A common technique for interpreting MWD/LWD and wireline data in 
environments with complicated geometry such as dipping beds is to employ a model which 
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computes estimates for the measurements as a function of the parameters of a hypothetical 
earth formation. Once model input parameters have been selected that result in a reasonable 
correlation between the measured data and the model data over a given depth interval, the 
model input parameters are then correlated to the actual formation parameters. This process 
is often referred to as "iterative forward modeling" or as "Curve Matching," and applying it 
in conjunction with the old assumptions, leads to errors because the volumes in which each 
measurement senses each variable have to be known in order to adjust the model parameters 
appropriately. 

[00107] The algorithms discussed in the previous sections can also be adapted for 
application to data acquired at non-zero relative dip angles. Selecting the background 
medium to be a sequence of layers having the appropriate relative dip angle is one method for 
so doing. 

TRANSFORMATIONS FOR A RESISTIVITY-DEPENDENT DIELECTRIC CONSTANT 
[00108] In the embodiments described above, both the dielectric constant and conductivity 
are treated as independent quantities and the intent is to estimate one parameter with minimal 
sensitivity to the other. As shown in Figure 1 , there is empirical evidence that the dielectric 
constant and the conductivity can be correlated. Such empirical relationships are widely used 
in MWD/LWD applications, and when they hold, one parameter can be estimated if the other 
parameter is known. 

[00109] This patent application shows that: 1) an attenuation measurement is sensitive to 
the conductivity in the same volume of an earth formation as the phase measurement is 
sensitive to the dielectric constant; and, 2) the attenuation measurement is sensitive to the 
dielectric constant in the same volume that the phase measurement is sensitive to the 
conductivity. A consequence of these relationships is that it is not generally possible to 
derive independent estimates of the conductivity from a phase and an attenuation 
measurement even if the dielectric constant is assumed to vary in a prescribed manner as a 
function of the conductivity. The phrase "not generally possible" is used above because 
independent estimates from each measurement can be still be made if the dielectric constant 
doesn't depend on the conductivity or if the conductivity and dielectric constant of earth 
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formation are practically the same at all points within the sensitive volumes of both 
measurements. Such conditions represent special cases which are not representative of 
conditions typically observed within earth formations. 

[00110] Even though two independent estimates of the conductivity are not generally 
possible from a single phase and a single attenuation measurement, it is still possible to 
derive two estimates of the conductivity from a phase and an attenuation measurement given 
a transformation to convert the dielectric constant into a variable that depends on the 
resistivity. For simplicity, consider a device such as that of Figure 3. Let the complex 
number w, denote an actual measurement (i.e., the ratio of the voltage at receiver 307 relative 
to the voltage at receiver 309, both voltages induced by current flowing through transmitter 
305). Let the complex number w denote the value of said measurement predicted by a 
model of the tool 320 in a prescribed earth formation 3 1 5 . For further simplicity, suppose the 
model is as described above in the section "REALIZATION OF THE 
TRANSFORMATION." Then, 

[0100] where the wave number k = k(cr, s{a)) = jicoju(a + icos{a)) , and the dependence of 
the dielectric constant s on the conductivity a is accounted for by the function e (<r) . 
Different functions s(a) can be selected for different types of rock. Let a P and a A denote 
two estimates of the conductivity based on a phase and an attenuation measurement and a 
model such as the above model. The estimates can be determined by solving the system of 
equations 

0=lw, \ -\w((T A ,£(CT p ))\ 

0 = phased) - phase(w(a P ,s{<7 A ))). 

[0101] The first equation involves the magnitude (a.k.a. the attenuation) of the measurement 
and the second equation involves the phase (a.k.a. the phase shift) of the measurement. Note 
that the dielectric constant of one equation is evaluated using the conductivity of the other 
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equation. This disclosed technique does not make use of the "old assumptions." Instead, 
the attenuation conductivity is evaluated using a dielectric value consistent with the phase 
conductivity and the phase conductivity is evaluated using a dielectric constant consistent 
with the attenuation conductivity. These conductivity estimates are not independent because 
the equations immediately above are coupled (i.e., both variables appear in both equations). 
The above described techniques represent a substantial improvement in estimating two 
resistivity values from a phase and an attenuation measurement given a priori information 
about the dependence of the dielectric constant on the conductivity. It can be shown that the 
sensitivity functions for the conductivity estimates a A and a P are S' and S" Respectively if 
the perturbation to the volume P 225 is consistent with the assumed dependence of the 
dielectric constant on the conductivity and a A = o> . 

[0102] It will be evident to those skilled in the art that a more complicated model can be used 
in place of the simplifying assumptions. Such a model may include finite antennas, metal or 
insulating mandrels, formation inhomogeneities and the like. In addition, other systems of 
equations could be defined such as ones involving the real and imaginary parts of the 
measurements and model values. As in previous sections of this disclosure, calibration 
factors and borehole corrections may be applied to the raw data. 

TRANSFORMATIONS FOR A RESISTIVITY-DEPENDENT DIELECTRIC CONSTANT 
USED IN CONJUNCTION WITH CONVENTIONAL PHASE RESISTIVITY VALUES 
[0103] The preceding section, "TRANSFORMATIONS FOR A RESISTIVITY- 
DEPENDENT DIELECTRIC CONSTANT," provides equations for o> and a A that do not 
make use of the "old assumptions." One complication that results from using said equations 
is that both a phase shift measurement and an attenuation measurement must be available in 
order to evaluate either a P or a A . In MWD/LWD applications, both phase shift and 
attenuation measurements are commonly recorded; however, the attenuation measurements 
are often not telemetered to the surface while drilling. Since bandwidth in the telemetry 
system associated with the logging tool is limited, situations arise where it is useful to have 
less precise resistivity measurements in favor of other data such as density, speed of sound or 
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directional data. The following reparameterization of the system of equations in the previous 
section accommodates this additional consideration: 

0=| w, \-\w(a A ,e(a p ))\ 

0 = phase(w x )-phase(w(a P , £(<?„))), 

[0104] where a P and a A respectively represent the reciprocals of the first and second 
resistivity values; the second equation involves a phase of the measured electrical signal and 
the first equation involves a magnitude of the measured electrical signal at a given frequency 
of excitation; a function *(•) represents a correlation between the dielectric constant and the 
resistivity and is evaluated in a P in the first and second equations; w, represents an actual 
measurement (e.g., the ratio of the voltage at receiver 307 relative to the voltage at receiver 
309, both voltages induced by current flowing through transmitter 305) in the form of a 
complex number, and a function w(<x, e ) represents a mathematical model which estimates 

w v More particularly, s(a P ) represents the transformation of dielectric constant into a 

variable that depends upon resistivity. 

[0105] The second equation evaluates phase conductivity a P and the dielectric constant 
correlation with the same phase conductivity which is consistent with sensing both resistivity 
and dielectric constant in substantially the same volume. The first equation evaluates 
attenuation conductivity a A and the dielectric constant correlation with the phase 
conductivity (not an attenuation conductivity) which is consistent with sensing the resistivity 
and dielectric constant in different volumes. The second equation allows a P to be 
determined by a phase shift measurement alone. The first equation is actually the same 
equation as in the previous section, but the value of a A that solves the equation may not be 
the same as in the previous section because the value of a P used in the correlation between 
the dielectric constant and the resistivity may not be the same. In other words, the phase 
conductivity a P can be determined without the attenuation measurement, but the attenuation 
conductivity a A is afunction of both attenuation and phase shift measurements. The valueof 
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a P resulting from these equations is the conventional phase conductivity cited in the prior art 
which is derived using the "old assumptions" that the phase shift measurement senses both 
the resistivity and dielectric constant within substantially the same volume. The value for a A 
is consistent with the fact the attenuation measurement senses the resistivity in substantially 
the same volume that the phase measurement senses the dielectric constant; however, the a A 
value is abetter approximation than the a A values in the previous section because the phase 
conductivity a P used to evaluate the correlation between the dielectric constant and the 
resistivity is derived using the "old assumptions." The above two equations are partially 
coupled in that the phase conductivity a P is needed from the second equation to determine 
the attenuation conductivity a A in the first equation. 

[0106] Since the phase shift measurement is typically less sensitive to the dielectric constant 
than the corresponding attenuation measurement, the equations given in this section for a P 
and a A provide a reasonable tradeoff between the following: 1) accurate resistivity 
measurement results, 2) consistent resistivity measurement results (reporting the same 
resistivity values from both recorded and telemetered data), and 3) minimizing amounts of 
' telemetered data while drilling. 
[0107] When the perturbation to the volume P 225 is consistent with the assumed 
dependence of the dielectric constant on the conductivity, the sensitivity functions for 
conductivity estimates a P and a A determined by the equations given in this section are: 

d ° AcrApAz 



I[S"]-co^-I[S'] 
da 



(57[5"]-57[5 r ]) 
^ da ) 
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[0108] where — represents the derivative of function s(a), and it was assumed for 
simplicity that <r A = a, . The sensitivity function for a P is the coefficient of AoApAz in the 
first equation, and the sensitivity function for a A is the coefficient of AcrA/?Az in the second 
equation. Note that the leading term S" in the sensitivity function of the first equation and 

the leading term in the sensitivity function of the second equation are approximately 
I[S'] 

the sensitivity function mentioned in the section, "TRANSFORMATIONS FOR A 
RESISTIVITY-DEPENDENT DIELECTRIC CONSTANT," and that the remaining terms 
are representative of the sensitivity error resulting from using the equations given in this 
section. In the second equation, the error term in the sensitivity for a A tends to be small 
becauseitisproportionalto (a>j^j . It can be shown that the error term for the sensitivity 
function of a a A estimate resulting from equations consistent with the old assumptions is 
proportional to and is thus a larger error than produced by the technique described in 



this section. As a result, the technique of this section represents a substantial improvement 
over the prior art. 

APPROXIMATIONS TO FACILITATE USE OF RAPIDLY-EVALUATED MODELS 
[01 09] Many of the techniques described in the previous sections are even more useful if they 
can be applied in conjunction with a simplified model of the measurement device. This is 
especially true for embodiments which use iterative numerical techniques to solve systems of 
nonlinear equations because the amount of computer time required to achieve a solution is 
reduced if details of the measurement device can be ignored. Such details include finite sized 
antennas 205, 207 and a metallic drill collar 203 . Models that do not include such details can 
often be evaluated rapidly in terms of algebraic functions whereas models including these 
details may require a numerical integration or other similar numerical operations. An 
alternative choice for solving such systems of equations is to store multidimensional lookup 
tables and perform inverse interpolation. An advantage of using the lookup tables is that the 
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model calculations are done in advance; so, results can be determined quickly once the table 
has been generated. However, 1) a large amount of memory may be required to store the 
lookup table; 2) it may be difficult to handle inputs that are outside the range of the table; 3) 
the tables may require regeneration if the equipment is modified, and 4) the tables themselves 
may be costly to generate and maintain. 

[01 10] The approach taken here is a compromise. One-dimensional lookup tables are used to 
renormalize each measurement so as to be approximately consistent with data from a 
simplified measurement device which has infinitesimal antennas and no metallic mandrel. 
The renormalized data are used in conjunction with a simpler, but more rapidly evaluated 
model to determine the conductivity and/or dielectric constant estimates. Higher dimensional 
lookup tables could be used to account for more variables, but this has proved unnecessary in 
practice. 

[0111] Suppose the function h 0 {a,e) represents a model that estimates the measurements as 
a function of the conductivity and dielectric constant which includes details of the tool that 
are to be normalized away (i.e. finite antennas and a metallic mandrel). Suppose h x {a, e) is a 
simplified model which is a function of the same formation parameters (i.e. the resistivity and 
the dielectric constant) but that assumes infinitesimal antennas and no metallic mandrel. The 
data shown in Figure 10 are for the medium spaced 2 MHz measurement described in 
conjunction with Figure 7. The first two columns of data resistivity (1/sigma) and dielectric 
constant (eps_rel) values are input into the respective models (which are not used in the 
algorithm, but are shown to illustrate how the table is generated). It was found that when the 
resistivity and dielectric constant are not both large, satisfactory results can be achieved by 
calculating the data as a function of the conductivity for only one dielectric constant, and in 
this embodiment, a relative dielectric constant of 35 was used. The third column (db_pt) 
contains the attenuation values for the simplified model evaluated as a function of 
conductivity, and the fourth column (dbjnan) contains the attenuation values for the model 
which explicitly accounts for the finite sized antennas and metallic mandrel. The fifth 
(degjpt) and sixth (deg_man) columns are similar to the third and fourth columns but for the 
phase shift instead of the attenuation. All data are calibrated to read zero if the electrical 
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parameters of the surrounding medium are substantially that of a vacuum (i.e., the air) 
(<r = 0,£ = 1 ). Specifically, columns 3 and 4 are: 

g' a = db_pt = 20 log 1 0(| \ (cr, e) I h { (0,1) |) 

ga = db_man = 201ogl0(| h Q (er,e)/ h 0 (0,l) |) 

[0112] Columns 5 and 6 are: 

g > =deg_pt = — arg^a, *)//*, (0,1)) 

p TV 
1 Rf) 

g =deg_man = — arg(A 0 (cr,£)/^ 0 (0,l)). 

[0113] In practice, calibrated measurements are used as values for db_man or phs_man (i.e. 
the ordinates in the one-dimensional interpolation). The corresponding values for db _pt or 
phs_pt result from the interpolation. Once the values for db_pt or phs_pt have been 
determined, relatively simple models such as 

( U Y &xy(ik 0 L 2 )(l-ik 0 L 2 ) 
1 ~{l 2 ) exp(^ 0 A)(l-^ 0 A) 

[01 14] can be parameterized as described in the several preceding sections and used to solve 
numerically for the desired parameters. The solutions can be obtained quickly, using robust 
and well known numerical means such as the nonlinear least squares technique. 
[0115] An alternative procedure would be to transform the values from the relatively simple 
model/ to corresponding db_man and deg_man values. This would produce equivalent 
results, but in conjunction with an iterative solution method, is clumsy because it requires the 
values from the simple model to be converted to db_man and deg_man values at each step in 
the iteration. 

EXTENDED APPROXIMATIONS TO FACILITATE USE OF RAPIDLY-EVALUATED 
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MODELS 

[0116] The embodiment described in above section, "APPROXIMATIONS TO 
FACILITATE USE OF RAPIDLY-EVALUATED MODELS," looses accuracy when the 
resistivity and dielectric constant are simultaneously large enough to cause the attenuation 
value to be less than its "air-hang" value. This region corresponds to values in Figure 8 with 
negative values on the horizontal axis. Embodiments of an extended technique addressing 
this shortcoming are described below. This extension is based on the observation that the 
difference between the mandrel and point-dipole values depends more strongly on the 
dielectric constant than the resistivity when the resistivity is large. 

[0117] When the attenuation and phase shift measurements are consistent with simultaneous 
large resistivity and dielectric constant values, the data in Figure 1 1 are used (instead of the 
data in Figure 10) to transform the measured values (deg_man, db_man) to the corresponding 
point-dipole values (deg_pt, dbjrt). Figure 12 is a plot of the difference between the point- 
dipole and mandrel values as a function of both the resistivity and dielectric constant. Two 
regimes are evident: 1) for low resistivity, the difference between the point-dipole and 
mandrel data depends more strongly on the resistivity than on the dielectric constant; and, 2) 
for high resistivity, the difference between the point-dipole and mandrel data depends more 
strongly on the dielectric constant than on the resistivity. 

[0118] In order to facilitate use of these extended approximations without operator 
intervention, point-dipole values consistent with both Figures 1 0 and 1 1 are computed. They 
are combined to form a weighted sum that biases the result heavily toward consistency with 
Figure 11 when appropriate and heavily toward consistency with Figure 10 when appropriate. 
Specifically, for the disclosed embodiment two terms are used in the weighted sum, 

g'a =C la f ia (g a ) + C 2a f 2a (g a ) 

. g' P =c {p f lp (g p ) + c 2p f 2p (g p ) 

[0119] where g' a and g' p represent the point-dipole values, g a and g p represent the 
measured (mandrel) values as defined under the heading, "APPROXIMATIONS TO 
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FACILITATE USE OF RAPIDLY-EVALUATED MODELS." The functions f u and f lp , 
defined numerically from data such as that of Figure 10, give corresponding point-dipole 
value as a function of the measured value for a first set of resistivity and dielectric constant 
values. The functions f 2a and f 2p , defined numerically from data such as that of Figure 11, 
give corresponding point-dipole value as a function of the measured value for a second set of 
resistivity and dielectric constant values. The coefficients c lfl , c 2a , c,„,and c 2p are chosen to 
produce values for g' a and g' p that approximate point-dipole values over the anticipated 
range of g a and g p . Satisfactory performance of the algorithm has been obtained with 

Ujgq -gqo) + eXp(~« | g a ~ gqO IMgoO - So) 

c = c, a =c lp =\-c 2a -\-c 2p - i + exp(-a|g o -g a0 l) 



[01201 where 



0, x < 0 
u(x) = \o.5,x = 0 
1 x>0 



a = 5 x 1 0" 5 /, / is the frequency in Hz 



[0121] and 



g a0 =J3H gp ) + r 



[0122] For each spacing and frequency the parameters fi and y can be determined 
empirically. For the medium 2 MHz measurement discussed above, fi = -0.047 dB and 
y = -0.159 dB are satisfactory. 

[0123] Figure 13 shows the results of applying this embodiment to the medium 2 MHz data. 
The dotted lines represent the measured (mandrel) data g a and g p . The dot-dash lines 
represent the corrected data g' m and g' p . Ideally, the corrected data would coincide exactly 
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with the actual point-dipole values which the solid lines represent. The difference between 
such ideal results and the results achieved using this embodiment are not significant formost 
practical applications. 

[01 24] It will be clear to those skilled in the art that further improvements in accuracy can be 
achieved if necessary. One way would be to use independent values for the coefficients c lfl , 

c 2a , c„ , and c 2p . Another alternative would be to add more terms to the weighted sums. If 

N terms are used, 

g' a =S c *fM 

N 

g' P =Z c i P Ms P ) 

[0125] where f ia and f ip represent the transformation between the mandrel and point-dipole 
data over prescribed trajectories in the g a ,g p plane. The embodiments described above in 
detail employs one trajectory determined by constant values of resistivity as in Figure 10 and 
another trajectory determined by constant dielectric constant values as in Figure 11. 
However, trajectories where neither variable is constant are also acceptable. 
MATHEMATICAL FORMULATION FOR ITERATIVE FORWARD MODELING USING 
THE FULL PARAMETER SPACE 

[0126] In the above section entitled, "ITERATIVE FORWARD MODELING AND 
DIPPING BEDS," dielectric constant and resistivity values are simultaneously determined 
using procedures known as "iterative forward modeling" or "curve matching." Those of 
ordinary skill in the art are generally familiar with a mathematical formulation for these 
procedures. A formulation is given here in order to define some notation that will be used in 
the next section. For simplicity, the electrical properties will be treated as piecewise constant 
functions of space. Consistent with the above text, the conductivity and dielectric values are 
treated as real-valued variables. This simplifying assumption results in dielectric losses being 
associated with the corresponding conductivity and in the association of conduction currents 
in quadrature with the electric field with the corresponding dielectric constant. 
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[0127] Let the vector 

w, = vv, +nv, = [w n ,...,w lw ] r 

[0128] denote a set of measured data where the superscript T denotes the matrix transpose. 
Each component of the vector w, is a complex number that represents a voltage ratio defined 
as explained in the descriptions of Figures 3 and 4. For convenience, it will be assumed that 
the number of measurements N = JK is the product of: 1) the number of depth values K 
being considered, and 2) the number of measurements / available at each depth. Suppose a 
model is available to estimate w, as a function of the vectors 
a =[cr,,...,cr i f 

[0129] and 

b=[b l ,...,b M f 

[0130] Each value cr,. and e t respectively represents a conductivity value and dielectric 
constant value within a prescribed region of space. Each value b, represents a spatial 
coordinate of a boundary separating two or more regions. Specific values of each vector 
<j,£,and b to be correlated with the parameters ofthe earth formation denoted by a\s \ 

and 6 * can be determined by solving the equation 

w, = F(a\ s',b*) 

[0131] where F represents a model used to estimate the measurement vector as a function of 
the vectors a ,£ , and b . The model F is itself a transformation (or function) that maps the 
vectors a.f.andMoa vector representative ofthe measurement. It is desirable for the 
Jacobian matrix of the transformation F to have a sufficient number of singular values 
within a domain near the point a\e\ and b* in order to ensure the problem is not 
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underdetermined. Many published numerical techniques for obtaining approximate solutions 
when a system of equations is fully, over, or under determined are known to those skilled in 
the art. Those of ordinary skill in the art are generally familiar with mathematical derivation 
and numerical implementation of mathematical models suitable for use as the vector function 
F. 

[0132] In some instances, it may be helpful to apply transformations to the measured and 
model data, and then solve the transformed system of equations. The formula 
) = h(f{g\ e' ,b* )) represents such transformed equations. The system of equations 
w, =F{U\e\b') is just the special case where the transformations G and H applied 
respectively to the measured and model data are both identity transformations. Examples 
where transformations other than the identity transformation are applied to the measured and 
model data are discussed below under the headings, "MATHEMATICAL FORMULATION 
FOR ITERATIVE FORWARD MODELING USING A REDUCED PARAMETER 
SPACE," and "USE OF PHASE AND/OR ATTENUATION RESISTIVITY DATA- 
MATHEMATICAL FORMULATION FOR ITERATIVE FORWARD MODELING USING 
A REDUCED PARAMETER SPACE 

[0133] In the previous section, the conductivity values are solved for along with the 
corresponding dielectric constants. Solving the system of equations 

[01341 may involve repeated evaluation of the Jacobian matrix for the transformation 
F which will be denoted by J(F) . This matrix contains the derivative of F with respect to 
each variable. After each evaluation, a singular value decomposition may be performed on 
J(F) . Evaluating and manipulating this matrix can be numerically intensive and time 
consuming. Use of the Cauchy-Reimann equations to relate derivatives of F i with respect to 
a j and £j helps reduce the cost of evaluating 1(F) , but does not reduce the dimensionality 
of this Nx(2L + M) complex-valued matrix where N represents the number of 
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measurements, L represents the number of conductivity and dielectric constant values and M 
represents the number of spatial coordinates being treated as variables. 
[0135] In many applications, the conductivity is the value of primary interest; consequently, 
the disclosed embodiment is a technique of determining values for ct* and b * without also 
determining e* , providing numerical efficiency. It should be evident that analogous 
procedures could be applied to estimate s' and b* without also determining a* . The steps 
of the disclosed embodiment are as follows: 

1. Apply a first transformation to the measured data such that the result, a 0 , is 
relatively sensitive to the conductivity and relatively insensitive to the dielectric 
constant: a 0 = Re(F 0 -1 (w, )) . 

2. Transform the unknown dielectric constant vector s which F depends on to a 
vector that is a function of the corresponding (but also unknown) conductivity 
vector a . In other words, s -» e(o) = [e, (cr, ), . . . , e L (a L )f . For convenience, a 
given dielectric value s i is assumed to be a function only of the corresponding 
conductivity value a, . Alternatively, some or all of the dielectric values could 
be assumed to be: 1) constants, independent of the conductivity; or 2) functions of 
any or all of the conductivity values associated with the vector a . 

3. Evaluate the model F(a ,s (a),b ) 

4. Apply the transformation in step 1 to F(a ,e {a),b ) 

5 . Identify values for a and b that cause the results of Steps 1 and 4 to agree, and 
correlate these values with a * and b * . 

[0136] A suitable choice for the transformation mentioned Step 1 is to use conductivity 
values determined as described above in "SINGLE MEASUREMENTS AT A SINGLE 
DEPTH," and described in more detail in, S. M. Haugland, "New Discovery with Important 
Implications for LWD Propagation Resistivity Processing and Interpretation," SPWLA 42 nd 
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Annual Logging Symposium, paper LL, June 17-20, 2001. In the embodiment disclosed 
here, Step 1 is realized by fitting each measurement (each component of the measurement 
vector w,), to a model of the measurement apparatus in a homogeneous medium and 
associating the conductivity of this medium with the result of the transformation. In terms of 
the above notation, 

w, =F 0 (Ia 0 +ine 0 ) 

[0137] where the NxN matrix J? is diagonal with each nonzero element being equal to the 
radian frequency of excitation of the corresponding measurement; 7 is the identity matrix; 
F 0 is the model of the measurement apparatus in the hypothetical homogeneous medium; the 
vector of length N , ct 0 , is the result of Step 1; and the vector of length N , £ 0 , is the 
dielectric constant vector associated with a 0 . No vector b was associated with F 0 because 
the medium is assumed to be homogeneous (however, an inhomogeneous background 
medium could be assumed without departing from the scope or spirit of this invention). In 
this embodiment, the system of equations w, = F 0 (Ia 0 + iQs 0 ) is uncoupled and can be 
expressed in scalar notation as follows: 

w lp =F Op (a 0p +ico p s 0p ); p = l,...,N 

[0138] The remaining steps are achieved by solving the system of equations: 
Re(F 0 -' (f{o* , e (a * ), b* ))) - v Q 

[0139] where Re(-) represents the real part of its argument. The inverse function F 0 ~ l has the 
properties: w, = F 0 (j? - )) and a 0 + iHe 0 = F 0 ~ ] (f o (Ia 0 + iHS 0 )). Realizations of the 
inverse function F 0 -1 applied to the model F and measured data w, maybe slightly different 
in order to account for differences between the measurement tool used to collect the 
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measured data and the properties of the measurement tool assumed in the model. 
Alternatively, the model F may be defined so as to directly account for these differences. 
INVASION EXAMPLE 

[0140] Instances where borehole fluids penetrate (invade) into the formation over time after 
drilling occur commonly. In some cases, the data acquired during the actual drilling process 
reflect invasion. In other cases, resistivity measurements from MAD (measurement after 
drilling) passes may reflect invasion even if the data acquired during drilling (the LWD data) 
did not. Invaded zones are often modeled by a piston-like cylinder surrounding the wellbore. 
Realizations of the vectors a, e,zndb to characterize the invaded region of the formation 
are:a = [<J X0 ,crJ,£= [s xo ,£, J , and b = [r t J where the subscripts xo and t respectively 
denote values associated with the invaded zone and virgin formation, and the radius of the 
invaded zone is . 

[0141] Given a set of measured data, one could solve directly for the two conductivity values, 
the two dielectric constants, and the invaded zone radius using the above mathematical 
formulation for iterative forward modeling using the full parameter space. 
[0142] In the above mathematical formulation for iterative forward modeling using a reduced 
parameter space, the measured data are transformed to an apparent conductivity vector cr 0 . 
The reduced sensitivity of a 0 to the dielectric constant is a motive transforming the unknown 
dielectric constant vector to a function of the unknown conductivity values: 

e s(a) = [e xo (<J xo ),£,(<? t )] ■ The remaining variables <r xo ,cj t and r ( are then estimated 
by selecting values for a and b such that computed estimates of the apparent conductivity 

values agree reasonably well with the apparent conductivity vector derived directly from the 

measurements. 

FREQUENCY DISPERSION 

[0143] The section entitled "MULTIPLE SENSORS AT MULTIPLE DEPTHS," mentions 
frequency dispersion is a complication that can arise if multiple excitation frequencies are 
used. One way to deal with frequency dispersion is to process data acquired using different 
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excitation frequencies separately. For example, if several excitation frequencies are used, 
then distinct measurement vectors would be formed for each excitation frequency, and 
distinct realizations of the vectors W X , and P would be computed for each exc.at.on 
frequency. Separate processing of multi-frequency data can be done in either the full or 
reduced parameter spaces. 

[0,44] Over the frequency range of interest in wireline induction and LWD propagation 
resistivity logging applications, frequency dispersion is more commonly observed on .he 
dielectric constant than the conduc.ivi.y. As a result, the technique for estimatmg the 
conductivity values in the reduced parameter space is expected to be relatively msens.t.ve.0 
dielectric frequency d.spersion even ifmul.i-ffequeneyda.a are simutaeously processed and 
no special measures are taken to account for (or suppress) dielectric dispersion effects. A 
variation of the embodiment disclosed above that is useful for further reducing the sens.fv.ry 

of the results to dielectnc frequency dispersion is .0 apply different realizations of funchon 

e (a) for each excitation frequency. 

,0145] A vanation of the above mathematical formulation for i.era,,ve forward modeling 
using the full pa— spaee that is useful to reduce the sensitivity of those results to 
dielectric fre,uencydispers,onis»oassocia.eadifferen.dielecmceons.an,vec.orw,,h each 

excitation frequency. Doing this results in the dielectric vector becoming the lx« mamx 



10146] where R represents the number of excitation frequencies. The unknown eonductivhy 
vecor could also be formulated in terms of a matrix to account for conductivity frequency 
dispersion. Such a matrix may be represented by: 
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,01471 in fact, mathematical formulation for iterative forward modeling using the full 
parameter spacing can be rewritten assumtng frequency dispersion on both variables by 
defining the variables in terms of appropriate matrices. In that case, the equatton 
w =F(o' X,b*) above becomes 

10148] and simafions where frequency dispers.on is no. accounted for can be handled by 
augmenting this system of equations with constraints causing some or all of the columns of 
either matrix W or T to be equal. Such constraint equations can be written in the form 
a = apr or Spq = e„. Defining unknown conductivity and d.electnc matrices m the 
JthemlalfoLaulfori^ 

account for frequency dispersion results in the system of equations: 
Re(F 0 -'(F(CT - ,f (o= - ),6')))=o- 0 

[0149] whichcanalsobeaugmentedbyconstraintsoftheforma,, = a„. These techniques 
reduce sensitivity to frequency dispersion and allow one to simultaneously process 
measurements from dtfferen, excttafion frequencies in the presence of frequency disperston. 
DROPPING THE PIECEWISE CONSTANT ELECTRICALP ARAMETER ASSUMPTION 
[01501 Useofthevector b as defined in the mathematical formmation for iterative modeling 
using the full parameter space and thereafter is an assumption that the electrical parameters 
are piecewise constant ftmcfions of space within the earth formation. Such simpl.fy.ng 
assumptions are commonly appl.ed in the art for convenience, however, more fafthful 
representations of spatial variations in the electrical properties are sometimes reqmred. For 
example, an abrupt transition between the virgin formation and invaded zone .s no. 
anticipated by mu.fiphase fluid flow models used to predict fluid profiles around the 
wellhore. Nonetheless, measured data areoftenfi.to electromagnetic modelswh.eh assume 
step boundaries between adjacent regions of the earth format.on. The pracfic. seems to be 
commonbeeause:l)electromagne.iemodelswhichexplic.,lyassumecontinuouslyvary.n g 
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(as a function of spatial variables) electrical parameters are not readily available; and 2) the 
spatial dependenceoftheresistivity and/or die.ectnc constant is typically no, known a pnon. 
It will be evident to those skil.ed in the art .ha, snch a simplifying assumption ,s optional. 
MAKING A DIELECTRIC ASSUMPTION WITHOUT APPLYING A 
TRANSFORMATION TO THE MEASURED DATA 

10151] in the mathematics formulation for iterative modeling using a reduced parameter 

space the dtelectric constants in the model F are assumed ,o be functions of the 

conductivities, and the remaining unknowns are estimated to give good agreement between 

the apparent conductivity values derived from ,he measured data and the apparen, 

conductivity values derived from data calculated using the mode, F For a measurement 

vector W, of length N , this method can be used to estimate no more than N conducnvny 

values. Upto 2N conductivity estima.es can be derived from the same measured data se, ,f 

,he apparent conductivity values are no. used and .he model is fi, directly te me measurement 

vector. Iftheassumeddielec.ricfunc.ionf 0 ) is a sufficiency accurate represent of 

theac.ualdielec.ric properties fitting *e model data directive .he measured datasubjecuo 
the dielectricassumptionwillgiveaccurate results. This embodimen.consis.sofso.v.ng.he 

system of equations 

w, =F(W',£ (<x*),&*) 

,01521 where .he consttain. equations mentioned above under .he headmg "FREQUENCY 
DISPERSION" can be applied to .he extent tha, frequency dispersion need nofbe accounted 
for. 

USE OF PHASE AND/OR ATTENUATION RESISTIVITY DATA 
,01531 Calculationofattenna«ionandphaseres,s.ivi.yvalues from individual measurements 
L dtscussed above under me headmgs, "TRANSFORMATIONS FOR A RESISTIVITY 
DEPENDENT DIELECTRIC CONSTANT" and "TRANSFORMATIONS FOR A 
RESISTIVITY DEPENDENT DIELECTRIC CONSTANT USED IN CONJUNCTION 
WITH CONVENTIONAL PHASE RESISTIVITY VALUES." Matching these reststivtty 
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values to similarly calculated resistivity values derived from model data is another 
embodiment that can be used. Suppose A. denotes a transformation that converts the 
measured voltage ratio value to a corresponding attenuation resistivity and P., denotesa 
transformation that converts measured voltage ratio value w„to the corresponding phase 
resistivity. Then, the real-valued vectors of length N , ? 0 (w, )and 2 0 (w, ) are respectively 
thephase and attenuation resistivity values. It is then possible to set up systems of equations 
of the form 



[0154] and 

A0>"i ) = a(f(?M (?•),*")) 

[01551 which can either be solved simultaneously or independently, depending on the 
situation. As discussed above, different versions of the transformations P 0 and A. can be 
usedontherightandkfthandsidesoftheaboveequatiousifdifTerenecsberiveentheaetual 
measurements and the measurement tool assumed in the mode, have not been accounted 
for by other means. It should be understood that any of the equations presented herem may 
be implemented inasoflware program executed by a processor of a compering system such 
as a computer at the surface or a "downhole" microprocessor. 

101561 The disclosed techniques provide accurate conductivity estimates that are 
independentofdieleetrieconstantproperties. Conductivity values serve several apphcattons 
such asWdetectingthepresence, absence or amount of a hydrocarbon, (ii)gnidmg a dn.lbtt 
within a productive zone, (iii) estimating poor pressure, and (iv) evaluating reservo.r and 
other geological features through correlation with logs and nearby wells. Two examples of 
applications for die.ectrie values are detecting the presence, absence or amount of a 
hydrocarbon and detecting vertical fractures. 
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,01571 The foregoing disciosvare and description of the various embodiments are 
inustrative and oratory thereof, and various changes in the descriptions, mo ^mg 
parameters and attributes of the system, the organization of the measurements, tr^sm r 
and receiver confrgnrations, and the order and timingofs.eps taken, as weU as m the detatis 

,00.58, AUbongb the present invention audits advantages have been descnbed mdeUd, t 
Lldbeunderstood.ha,variouschanges,substi,utionsandmte m a,ionscanbemadehercm 

claims. 



